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1  INTRODUCTION 


This  report  describes  an  approach  for  calculating  the  bistatic  microwave  reflectivity  of  the  sea 
surface.  This  report  will  only  address  low  grazing  angles,  as  encountered  with  shipboard  radar  systems, 
but  include  both  in-plane  and  out-of-plane  geometries.  Higher  grazing  angles  as  well  as  airborne  or  space- 
based  radars  will  need  additional  models. 

In  the  forward  scatter  region  the  model  is  based  on  the  work  of  Beckman  and  Spizzichino  [1] 
supplemented  with  later  work  on  shadowing  [2]  and  polarization  effects  [3].  Outside  the  forward  scatter 
glistening  zone,  the  model  is  heuristic  and  based  on  similarities  between  measured  monostatic  and  bistatic 
reflectivities.  The  main  objective  of  the  model  described  here  is  to  evaluate  the  magnitude  of  jamming 
reflections  into  the  radar  main  beam  and  sidelobes  from  the  sea  surface.  Such  jamming  reflections, 
sometimes  referred  to  as  “hot  clutter”,  can  degrade  the  performance  of  adaptive  jamming  cancellation 
systems  such  as  side  lobe  cancellers.  The  model  would  also  apply  to  the  evaluation  of  the  performance  of 
bistatic  radar  systems. 

Under  standard  propagation  conditions,  bistatic  reflections  via  the  sea  surface  for  a  remote 
transmitted  signal  into  a  receive  antenna  will  exist  at  all  azimuth  angles  out  to  the  4/3  earth  radar  horizon 
associated  with  either  transmitter  or  receiver  heights.  During  ducting  conditions,  reflections  can  occur 
from  longer  ranges  dependent  on  the  refractive  properties  of  the  atmosphere.  The  effects  of  ducting  could 
most  likely  be  accounted  for  in  the  model  using  a  procedure  similar  to  that  proposed  by  Dockery  [4]  for 
monostatic  sea  clutter.  This  approach  relies  on  propagation  predictions  obtained  from  a  program  such  as 
TEMPER  [5].  This  generalization  is  not  included  in  the  present  work. 

2  BISTATIC  RECEIVED  POWER 

The  bistatic  power  reflected  from  a  given  radar  resolution  cell  is  proportional  to  the  product  of  the 
resolution  cell  area  (Ac)  and  the  normalized  bistatic  reflectivity  (cr^)  of  the  sea  surface  for  a  given 

geometry.  The  bistatic  reflectivity  has  the  non-dimensional  unit  of  bistatic  radar  cross-section  per  unit 
surface  area  (mVm^).  The  received  bistatic  power  further  takes  into  account  the  transmitted  power,  the 
gain  of  both  transmit  and  receive  antennas,  and  the  one-way  propagation  factors  from  the  transmit 
antenna  to  the  reflecting  area  on  the  surface  and  from  here  to  the  receive  antenna. 

For  non-cooperative  scenarios,  such  as  noise  jamming,  the  resolution  cell  area  is  defined  by  the 
intersection  of  transmit  and  receive  beams.  Pulse  jamming  will  additionally  need  to  consider  the  spatial 
extent  of  the  pulse  length.  For  cooperative  systems,  such  as  bistatic  radar,  the  resolution  cell  area  is 
further  limited  by  the  surface  area  contained  within  iso-range  and  iso-doppler  contours  associated  with  the 
range  and  Doppler  resolution  of  the  system.  In  cases  where  the  resolution  cell  area  is  so  large  that  the 
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reflection  coefficient  varies  spatially  within  the  cell,  the  cell  must  he  subdivided  and  contrihutions  non- 
coherently  summed  to  calculate  the  total  average  power. 


The  superposition  of  three  physical  processes  descrihes  histatic  sea  reflectivity.  These  are  (1) 
coherent  scatter  (reflection)  from  the  specular  region,  (2)  strong  diffuse  scatter  from  the  glistening  region, 
and  (3)  wide  angle  diffuse  scatter  from  the  sea  surface  outside  the  glistening  region.  Figure  1  shows  these 
three  regions  for  an  antenna  located  on  a  Navy  ship.  The  results  apply  only  to  the  sea  surface  inside  the 
4/3  earth  horizon  for  both  transmitter  and  receiver.  Figure  la  shows  the  volume  of  the  three  physical 
processes  described  earlier.  Figure  lb  shows  an  example  of  the  glistening  zone. 


Figure  la:  The  three  physical  processes  for  histatic  sea  reflectivity 
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^  Glistening  Zone  ~  4  deg  W  j 

Figure  lb:  An  example  of  a  glistening  zone 


Specular  forward  scatter  from  the  sea-surface  will  result  in  coherent  multipath  interference  while 
diffuse  forward  scatter  will  produce  non-coherent  scatter.  This  report  discusses  each  of  these  three  forms 
of  scatter  and  combines  them  into  a  MATLAB®  script  used  to  calculate  the  surface  reflectivity  for 
arbitrary  geometries.  Each  of  the  three  cases  uses  a  different  mathematical  model. 

2.1  Specular  Scatter 

Specular  scatter  comes  from  a  small  region  on  the  surface  around  the  specular  point  (see  Fig.  la)  at 
which  the  incidence  and  reflection  angles  are  equal  and  located  in  the  vertical  plane  including  both 
transmitter  and  receiver.  The  effect  of  standard  refraction  by  the  atmosphere  is  taken  into  account  by 
using  an  effective  earth  radius  Re  -  8500  km  instead  of  the  actual  radius  of  the  earth.  The  range  loss 
associated  with  the  total  path  length,  from  the  source  to  the  specular  point,  and  from  here  to  the  receiver, 
together  with  the  coherent  reflection  coefficient  of  the  surface  determines  the  complex  amplitude  of  the 
received  specular  reflected  signal.  In  addition,  the  gain  of  the  transmitter  and  receiver  antennas  must  also 
be  taken  into  account. 


The  coherent  reflection  coefficient,  pc  is 

Pc  =Po-Ps  D 
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In  Eq.  (1),  is  the  Fresnel  reflection  coefficient  for  sea  water  for  either  horizontal  or  vertical 
polarization,  is  the  additional  rough-surface  specular  reflection  coefficient,  and  D  is  the  divergence 
factor  (see  p.245  in  [1]).  For  the  histatic  model  addressed  here,  the  divergence  is  assumed  to  he  unity. 

The  Fresnel  reflection  coefficients  can  he  calculated  using  classical  EM  (electromagnetic)  equations 
(see  p.  218-219  in  [1]).  The  Fresnel  reflection  coefficient  for  horizontal  polarization  is: 


H 


Po 


sin  Y  -  -  cos^  y 

sin  y  +  ^/y^'^-cos^ 


and  for  the  vertical  polarization. 

sin  y  -  Jy^  -cos^  y 

Po  ~  ^  I 

Y  sin  y  +  yJY  -  cos  y 


(2) 


(3) 


In  these  equations,  Y  is  the  electrical  impedance  of  the  surface: 


Y  = 


(4) 


where  //„  is  the  relative  permeability  ( »  1 )  and  Src  is  the  complex  dielectric  constant  defined  as: 


^rc  =^r-  JSOAO- 


(3) 


Here  Sy  is  the  relative  permittivity  of  the  sea  surface  assumed  to  he  80  (see  p.219  in  [1]),  A.  is  the 
wavelength  in  meters  and,  cris  the  conductivity  of  the  surface  in  mho/m. 

The  angle  y  is  the  grazing  angle  of  the  incident  field.  For  specular  scatter,  the  grazing  angle  for  the 
incident  field  is  equal  to  the  grazing  angle  of  the  reflected  field.  As  the  sea  surface  gets  rougher,  the 
specular  component  reflects  smaller  amounts  of  the  EM  energy.  The  rough  surface  specular  reflection 
coefficient,  ps  quantifies  this  phenomenon.  The  rough  surface  specular  reflection  coefficient  depends  on 
the  root-mean-square  (rms)  wave  height,  ct/,,  the  wavelength  of  operation.  A,,  and  the  grazing  angle,  y.  For 
a  Gaussian  sea  surface  model,  the  rms  rough  surface  specular  reflection  coefficient  is  given  hy  (eq.246  in 
[1]): 


Ps  =  exp 


4;Tcr^  sin  y 

.  I 


(6) 
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Using  Eq.  (6),  Fig.  2  shows  the  specular  coefficient  for  various  sea  states.  When  sea  state  is  zero,  the 
reflection  coefficient  is  one  and  it  drops  rapidly  as  the  sea  state  increases.  This  plot  is  at  X-hand  (10 
GHz).  Section  2.2. 1  describes  the  Gaussian  sea  surface  model  in  more  detail. 


Figure  2:  Specular  reflection  coefficient  vs.  grazing  angle  for  different  sea  states 
2.2  Diffuse  Scatter 

Diffuse  scattering  occurs  over  an  extended  region  of  the  sea  surface.  The  primary  contrihution  is 
from  a  region  known  as  the  “glistening  surface”  (p.250  in  [1]).  The  glistening  surface  is  assumed  to 
consist  of  small  wave  facets  (presumed  larger  than  radar  wavelength),  which  act  as  mirrors  and  each  of 
which  is  oriented  such  that  the  transmitter  reflects  an  incoming  EM  field  towards  the  receiver.  Figure  3 
illustrates  such  diffuse  scattering.  In  Fig.  3,  any  wave  facet  with  the  correct  slope  and  orientation  reflects 
the  EM  energy  toward  the  receiver.  If  the  slope  is  not  correct,  then  this  energy  is  scattered  in  other 
directions. 
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Transmitter 


Figure  3:  Diffuse  scattering 

2.2.1  Sea  Surface  Model 

The  model  assumes  that  the  rough  sea  surface  lies  in  the  x-y  plane  with  the  height  of  the  rough 
surface  in  the  z-direction.  The  rough  sea  surface  is  then  modeled  as  a  two-dimensional  Gaussian  random 
surface  according  to  Beckmann  and  Spizzichino  (p.  80,  eq.(2)  in  [1]).  The  term  defines  the  sea 

surface  where  is  the  wave  height  relative  to  the  mean  at  coordinate  (x,  _y) .  The  resulting  Normal 
prohahility  density  function  of  the  wave  height  is  defined  hy  the  rms  (root-mean-square)  value  of  the 
wave -height,  ct/,,  and  is  given  hy 


p(C) 


1 


f 

exp 

V 


2cT;?J 


C7J 


The  rms  wave  height  ct/,  is  14  of  the  significant  wave  height  which  is  defined  as  the  average  of  the 
one-third  highest  waves  at  the  given  sea  state.  Table  1  shows  the  definition  of  sea  states  as  provided  hy 
the  World  Meteorological  Organization  (WMO)  and  the  related  significant  wave  heights,  rms  wave 
height,  and  correlation  distance. 

The  autocorrelation  function  for  the  two-dimensional  Gaussian  random  surface  is  assumed  to  he 
Gaussian  and  is  defined  hy  (p.81  in  [1]): 
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CCz")  =  exp 

V 


(8) 


In  eq.  (8),  T  represents  the  “correlation  distance”  of  the  waves.  For  this  model,  the  value  of  T  is  usually 
much  larger  than  the  wavelength  of  the  radar.  The  model  assumes  that  the  correlation  distance  is  the 
same  in  all  directions.  Table  1  also  includes  the  correlation  distance  T  assumed  for  the  different  sea  states. 
An  example  of  the  Gaussian  sea  surface  model  based  on  equations  (7-8)  is  shown  in  Fig.  4. 


Figure  4:  Simulated  Gaussian  sea  surface  model 


2.2.2  Diffuse  Scatter  Reflectivity 

As  previously  described,  the  direction  and  amplitude  of  the  reflected  wave  assumes  that  the  sea 
surface  between  the  receiver  and  transmitter  is  entirely  composed  of  small  plane  elements  (also  referred 
to  as  facets).  Each  of  these  facets  acts  as  a  mirror.  The  distribution  of  the  slopes  of  these  mirrors  will  have 
a  normal  distribution  with  a  zero  mean  and  a  standard  deviation  defined  as  the  rms  wave  slope,  go-  Based 
on  the  Gaussian  surface  model  of  the  sea,  it  can  be  shown  that  the  rms  wave  slope  is  given  by  the  rms 
value  of  the  wave  height  ct/,  and  the  correlation  distance  T  according  to  (p.251  in  [1]): 


^0  =tan/?o  = 


T 


(9) 


This  equation  also  defines  the  angle,  /?o,  corresponding  to  the  rms  slope.  Due  to  the  symmetry  of  the 
sea  surface  model  the  rms  wave  slope  is  the  same  in  all  directions.  The  “glistening  surface”  is  the  region 
from  which  the  transmitter  reflects  significant  power  of  the  incident  EM  wave  into  the  receiver  by  facets 
with  the  appropriate  slope  and  direction.  Table  1  lists  the  rms  wave  slopes  for  the  different  sea  states. 
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Table  1:  Sea  States  Wave  Height,  Correlation  Distance,  and  RMS  Wave  Slope 


WMO 

Sea 

State 

Significant 

Wave 

Height(m) 

Characteristics 

RMS  Wave 

Height 

Oh  (m) 

Correlation 

Distance 

T(m) 

RMS  Wave  Slope, 

? 

0 

<0.04 

Calm 

-0.01 

0.40 

0.05 

1 

0.04-0.1 

Calm  (ripples) 

0.03 

0.47 

0.12 

2 

0. 1-0.5 

Smooth 

0.11 

1.61 

0.14 

3 

0.5-1.25 

Smooth 

0.29 

3.92 

0.15 

4 

1.25-2.5 

Moderate 

0.59 

7.41 

0.16 

5 

2.5-4.0 

Rough 

1.03 

11.42 

0.18 

6 

4.0-6.0 

Very  Rough 

1.61 

14.68 

0.22 

7 

6.0-9.0 

High 

2.37 

18.96 

0.25 

2.2.3  Normalized  Bistatic  RCS 

In  Beckmann  and  Spizzichino,  [1]  an  equation  was  derived  for  the  bistatic  RCS  based  on  an 
evaluation  of  the  statistical  properties  of  the  wave  slope  for  a  given  transmitter/receiver  geometry. 
However,  this  equation  did  not  consider  the  effects  of  surface  reflectivity,  polarization,  and  shadowing. 
To  include  all  of  these  effects,  the  product  of  the  three  terms  defines  the  normalized  bistatic  RCS  as: 

<  y)  =  G(P{x,  y))  ■  S  (x,  y)  ■  \p{x,  yf 

In  Eq.  (10),  G{j8(x,  y))  accounts  for  the  wave-slope  statistics  (the  classical  Beckmann  and  Spizzichino 

equation),  S(x,  y)  accounts  for  shadowing,  and  p(x,  y)  represents  the  reflectivity  and  polarization  effects  of 
the  individual  wave  facets. 

The  term  G{/3(x,  y))  accounts  for  the  wave-slope  statistics  at  a  surface  cell  whose  location  is  given 
by  its  X-  and  y-  coordinates.  Assuming  that  the  rough  sea  surface  is  modeled  as  a  two-dimensional 


Gaussian  random  surface,  Beckmann  and  Spizzichino  have  obtained  the  expected  value  of  the  reflectivity 
as  a  function  of  the  required  slope  of  the  facet,  P ,  and  the  rms  wave-slope,  Po ,  as  (p.252  in  [1]): 


G{P) 


1 


tan^  {P) 
tan"(/?o) 


(11) 


The  required  slope  P ,  is  the  angle  between  the  bisector  angle  of  the  incident  and  scattered  rays  and  the  z- 
axis  of  the  surface  as  shown  in  Figure  5. 


Transmitter 


To  calculate  the  bisector  angle,  /?,  assume  that  the  receiver  is  located  at  the  origin  with  a  finite  height 
of  ha  and  the  transmitter  is  located  at  a  distance  d  away  at  height  hj.  The  x-y-z  coordinates  are  defined 
such  that  the  reflecting  cell  and  the  receiver  are  always  located  on  the  x-axis.  The  angle  Yr  and  yj  are  the 
grazing  angles  of  the  receiving  and  transmitting  paths  with  the  x-y  plane.  The  bisector  angle  can  be 
determined  in  terms  of  the  receiver  grazing  angle,  the  transmitter  grazing  angle  and  the  sum  of  the  off 
axis  angle  of  the  transmitter  02  and  off  axis  angle  of  the  receiver  0i  as 


JcoP  2cos(xr  )cos(^^  ) 

tan  P  =  - 


cos 


ie)- 


cos 


iru) 


sm 


(/«)  +  sin(Xr) 


(12) 
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where  the  angle  d  is  the  sum  of  the  angles  0i  and  62- 

In  the  case  of  a  curved  earth,  the  grazing  angles  are  calculated  with  respect  to  the  horizontal  at  the 
reflecting  cell  as  discussed  hy  Blake  (pp.254-257  in  [6]).  Figures  6a,  h  and  c  show  the  curved  earth 
geometry  for  the  three  paths  -  d  (surface  path  between  receiver  and  transmitter),  r/j- (surface  path  between 
transmitter  and  reflecting  cell)  and  (surface  path  between  the  receiver  and  the  reflecting  cell). 

Considering  the  geometry  of  Figure  6a,  the  angle  84  relates  the  height  of  the  transmitter  and  the 
height  of  the  receiver  by 


(Re  +/t«)xsin(6>^  +^) 

/v  = - — 


R 


(13) 


where  the  angle  (f)  is  measured  at  the  center  of  the  earth  and  is  equal  to  —  . 

The  path  lengths  Ri  and  R2  together  form  the  indirect  path  while  the  path  length  R4  is  the  length  of  the 
direct  path.  The  path  lengths  Rj  and  R2  can  be  calculated  using  the  geometry  shown  in  Figures  6b  and  6c 
as 


R.= 

R2  = 


Ih^  +4xR^x(^R^  +hj.)  X  sin^ 


[2) 


where 

Re 

K 


(14) 


(15) 


An  indirect  path  between  a  receiver,  transmitter  and  a  reflecting  cell  will  be  possible  only  by 
establishing  a  line  of  sight  between  the  three  points. 

Finally,  the  grazing  angles  Yr  and  yj  can  be  calculated  using  the  geometries  in  Figures  2b  and  c 
respectively.  Since  the  distances  in  the  y-direction  are  small,  the  angle  9  can  be  accurately  determined 
using  flat  earth  geometry. 
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Figure  6a:  Curved  Earth  Geometry  when  the  Surface  Path  (d)  between  Receiver  and  Transmitter  is 

considered 


Receiver 


Figure  6h;  Curved  Earth  Geometry  when  the  Surface  Path  (dR)  between  Receiver  and  Reflecting  Cell  is 

considered 
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iTransmitter 


Reflecting 

Surface 


Figure  6c:  Curved  Earth  Geometry  when  Surface  Path  {dr)  between  Transmitter  and  Reflecting  Cell  is 

considered 


The  term  /^(x,  y)  calculates  the  reflectivity  of  the  wave  facet  (reflecting  cell)  located  at  coordinates 

X,  y.  This  term  includes  the  change  in  polarization  due  to  reflection  of  a  tilted  surface.  The  reflectivity  of 
a  wave  facet  is  given  hy  (eq.  18  in  [3]): 


Pw 

Phv  = 

PVH  = 

Phh  = 


-a  + 


V 


-Po 


u 


cos(;i'g)cos(/j.)sin^  {6) 
sin^  {2a)  ^ 

aj  cos(;K^)sin(0) 


H 


^  n 

- a 

u  . 


y 


sin^  (2a) 


H 


H 


-Po 


K 


\2 


-a  + 


cos 


a  + 
) 


-  y  sin^(2a) 
cos(f^)cos(fj.)sin^  (^) 
sin^  (2a) 


Po 

f 


n 

U 

f 


\ 


-a 


n 

U' 


-a 


n 


-a 


(16) 
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where  the  first  subscript  is  the  polarization  of  the  scattered  signal  and  the  second  subscript  is  the 
polarization  of  the  incident  signal  and  the  parameters  and  are  the  Fresnel  reflection  coefficients  as 

described  in  Section  2.1  .In  section  2. 1  and  were  calculated  at  the  grazing  angle,  y,  which  was 
the  grazing  angle  with  respect  to  the  horizontal.  In  this  case,  the  grazing  angle  will  be  determined  from 
the  angle  a  which  is  the  angle  between  the  bisector  and  either  the  incident  or  the  reflected  ray.  It  can  be 
calculated  in  terms  of  the  grazing  angles  and  the  off-axis  angle  of  the  transmitter  as: 


cos 


(a)  = 


/l-cos(f^)cos(fy  )  008(6*) -f  sin  sin  (fj.) 


(17) 


Note,  that  Eq.  (16)  considers  a  general  case  where  polarization  of  the  transmitted  beam  can  change  as 
it  scans  in  azimuth  (electronically  scanned  array).  For  example,  at  boresight,  a  vertical  beam  is  truly 
vertical  but  as  the  main  beam  scans,  some  of  the  energy  will  leak  into  the  orthogonal  polarization 
(horizontal).  Thus,  the  second  terms  in  Eq.  (16)  are  to  include  the  effects  of  these  cross  pol  terms.  The 
geometrical  parameters  a2  and  as  are: 

02  =  cos(^^ )  sin(^2- )  +  sin(;'^)cos(^2)‘^o®^^) 

03  =  cos(;K7-)sin(;Kg) -I- sin(;K2’)cos(;K^)cos(0) 


2.2.4  Shadowing  Effects 

Neglecting  shadowing  effects,  when  considering  scattering  of  EM  waves  from  rough  surfaces,  can 
result  in  errors,  as  the  power  due  to  diffuse  scattering  will  be  overestimated.  Several  papers  such  as  Smith 
[7]  and  Wagner  [2]  have  addressed  this  issue.  According  to  these  authors  the  simplest  way  of  including 
shadowing  of  a  rough  surface  is  to  introduce  a  shadowing  factor  S(x,y)  at  the  location  of  the  reflecting 
facet. 

For  there  to  be  no  shadowing,  a  given  facet  has  to  be  “visible”  to  both  the  receiver  and  the 
transmitter.  Otherwise,  the  reflecting  cell  will  be  shadowed  and  will  not  contribute  to  diffuse  scattered 
power.  Figure  7  shows  a  ray  being  shadowed.  This  discussion  implies  the  assumption  that  geometrical 
optics  provides  a  valid  mathematical  approach.  Presently,  there  is  no  known  approach  that  takes 
diffraction  into  account.  The  equations  to  calculate  shadowing  have  been  derived  from  Smith  (eqs.20-25 
[7])  and  Bourlier  (eqs.5a,49  in  [8])  and  are  one-dimensional  in  nature.  These  equations  are  dependent 
only  on  the  slopes  of  the  rays  from  the  receiver  and  the  transmitter  but  are  independent  of  the  heights  of 
the  waves: 
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(19) 


S{x,y) 


(1  +  etfjvj^ ))  X  (1  +  erfiy^ )) 
4^1  +  +  Ay-  ^ 


where 

_  |tan(r, 

*  ^/2  tan  Pq 

&'’"-v^4^erfc{v^) 


In  this  equation,  k  corresponds  to  either  the  transmit  path  (7)  or  the  receive  path  {R). 


Receiver 

Rx 


Transmitter 

Tx 


Figure  7:  Effect  of  Shadowing 


Using  Eq.  (19-20),  the  shadowing  factor  is  plotted  vs.  incident  grazing  angle  for  three  different 
transmit  grazing  angles  of  1°,  5°  and  10°.  These  plots  show  that  shadowing  decreases  with  increasing 
transmit  and  receive  grazing  angles.  This  means  that  there  is  higher  prohahility  of  a  given  resolution  cell 
being  “seen”  hy  both  the  transmitter  and  receiver.  Also  at  higher  sea  states,  the  shadow  factor  is  larger. 
Once  again,  higher  sea  states  mean  increased  wave  heights  and  hence  an  increased  probability  of  a 
surface  cell  being  shadowed. 
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Transmit  Ray  Grazing  Angle  =  1°  Transmit  Ray  Grazing  Angie  =  5° 


Figure  8:  Shadowing  function  vs.  incident  grazing  angle  for  three  different  transmit  ray-grazing  angles 
2.2.5  Wide  Angle  Scatter 

In  addition  to  the  forward  reflection  caused  by  specular  reflection  and  scattering  from  the  glistening 
surface,  diffuse  scattering  occurs  at  wide  angles  outside  the  glistening  surface.  This  is  the  result  of, 
capillary  waves,  spray,  and  breaking  waves,  and  the  mechanisms  are  similar  to  those  causing  monostatic 
sea  clutter  returns.  In  this  bistatic  model,  the  wide-angle  scattering  is  determined  simply  by  calculating 
the  monostatic  RCS  for  the  smaller  of  the  two  grazing  angles  on  transmit  and  receive  -  yr  or  The 
maximum  of  the  contributions  from  the  glistening  surface  or  the  wide-angle  scattering  determines  the 
normalized  bistatic  RCS  for  any  azimuth  angle.  The  normalized  bistatic  RCS  for  wide-angle  scattering  is 
calculated  using  an  empirical  equation  [9]  based  on  the  Nathanson  tables  of  experimental  data  [10].  The 
bistatic  RCS  for  the  vertical  or  horizontal  polarizations  for  the  given  reflecting  cell  is: 
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(Jy  fj  (dB)  =  q  +  Cj  X  logjo  (sin  [min  )])  + 


(27.5-I-C3  X 

Lmin(^^,^^)j 

)x  logic/ 

1  +  0.95X 

Lmin(^^,^^) 

C4  X  (l  +  55)2+0.085x[min(r«,r7-)]+0.033xS5  +  C5  •  [min  ,  /j  )] 


(21) 


where  /t  and  /r  are  the  grazing  angles  (in  deg)  determined  in  Section  2.2.3  ,  SS  is  the  sea  state  and /is 
the  operating  frequency  of  the  radar  in  GHz. 

Then  for  wide  angle  scatter 


=  maxi  cr„ 


(cry,cr  / 


(22) 


The  five  parameters  for  the  calculation  of  the  RCS  for  the  horizontal  and  vertical  polarizations  are: 

Table  2:  Constants  Used  in  Empirical  Sea  Clutter  Model 


Constants 

Polarization 

Horizontal 

Vertical 

Cl 

-73.0 

-50.79 

C2 

20.78 

25.93 

C3 

7.351 

0.7093 

C4 

25.65 

21.58 

Cs 

0.0054 

0.00211 

Figure  9  plots  the  empirical  sea  clutter  model  equations  and  compares  them  with  the  points  from 
Nathanson’s  [10]  experimental  data.  The  results  shown  are  for  S-hand  (3.0  GHz)  for  sea  states  ranging 
from  zero  to  six  and  grazing  angles  varying  from  0. 1  -90  deg.  Analysis  has  shown  that  the  mean  absolute 
deviation  between  Nathanson’s  tables  and  Eq.  (21)  is  less  than  2.3  dB  [11]. 
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Figure  9:  Comparison  of  empirical  vs.  experimental  data  for  (a)  horizontal  and  (b)  vertical  polarization 
2.2.6  Examples 

Using  the  equations  summarized  in  the  earlier  sections,  some  numerical  examples  are  shown  in  the 


following.  These  results  were  obtained  using  the  MATLAB®  program  described  in  the  Appendix.  In 
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addition  to  the  total  bistatic  surface  reflectivity,  results  are  also  obtained  showing  the  time  delay 
distribution  of  power  reflected  into  the  receiving  antenna.  Such  time  delays  play  an  important  role  in 
determining  the  maximum  jamming  cancellation  that  can  be  achieved  with  a  given  adaptive  architecture. 

Figure  10  plots  the  bistatic  reflectivity  of  the  glistening  zone  over  a  surface  area  of  10  km  X  300  m. 
The  receiver  antenna  is  at  a  height  of  20  m  while  the  transmitter  has  a  height  of  225  m  above  the  surface. 
The  transmitter  has  an  isotropic  gain.  The  receiver  has  a  beamwidth  of  2°  X  2°  at  S-band  (3  GFlz).  All 
results  are  normalized  to  the  gain  of  the  receiver.  The  sidelobes  of  the  receive  antenna  are  40  dB  below 
the  main  beam.  The  sea  state  is  three.  In  the  example  shown  in  Fig.  10,  neither  shadowing  nor  wide-angle 
scatter  is  included. 


Reflectivity  of  Surface  (dB) 


Y-direction  (m)  MinPiotvaiue 


(OB) 

Figure  10:  Example  of  bistatic  reflectivity  from  glistening  zone 

Since  the  glistening  area  is  quite  large,  reflections  from  the  surface  have  a  distribution  in  time.  Power 
from  surface  can  reach  the  receiver  at  from  times  as  short  as  a  couple  of  nanoseconds  to  times  in  hundreds 
of  nanoseconds  and  even  microseconds.  Figure  11a  shows  an  example  of  the  distribution  of  time  over  the 
surface  by  drawing  isolines  (in  time)  of  10,  20,  30  and  40  ns.  In  this  example,  the  reflections  from  the 
edge  of  the  glistening  zone  will  take  more  than  40  ns  to  reach  the  receiver.  To  draw  a  plot  of  power  vs. 

time,  sum  the  reflected  power  bound  by  regions  of  equal  time.  Using  this  technique.  Fig.  1  lb  plots  of  the 
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power  distribution  vs.  time  and  cumulative  power  distribution  vs  time.  A  similar  technique,  using 
isolines  in  range,  allow  forming  a  plot  power  distribution  vs.  range.  Appendix  B  shows  an  example  of  a 
range  plot. 
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Figure  11:  (a)  Isolines  to  indicate  regions  of  equal  time  (b)  Power  distribution  vs.  time 
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In  this  next  example  (Fig.  12a  and  12b),  the  wide-angle  scatter  is  also  included.  The  surface  is 
extended  to  +10  km  in  the  y-direction  to  include  the  entire  wide-angle  scatter.  With  the  wide-angle 
scatter  included,  reflected  power  returns  arrive  with  delays  as  long  as  10  ms!  The  returns  with  shorter 
delays  are  from  the  glistening  region  while  the  returns  with  longer  delays  are  due  to  the  wide-angle  scatter 
(Fig.  12  b) 


(a) 
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Figure  12  :  (a)  Bistatic  surface  reflectivity  plot  (b)  Time  vs.  power  plots 


Finally,  the  effect  of  shadowing  is  shown  via  the  surface  reflectivity  plots  in  Fig.  13.  In  this  case,  the 
focus  is  only  at  the  glistening  region. 


Figure  13:  Comparing  the  change  in  the  surface  reflectivity  in  the  glistening  region  with  and  without 

shadowing 
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2.2. 7  Normalized  Radar  Cross-Section:  Ducting  Conditions 

Given  the  complexity  of  modeling  bistatic  sea  clutter  reflectivity,  as  described  in  the  previous 
Sections,  with  its  superposition  of  three  different  physical  processes,  accounting  for  the  effects  of  ducting 
could  be  complicated.  Attempting  to  calculate  bistatic  sea  reflectivity  under  ducting  conditions  in  a 
manner  analogous  to  that  used  for  monostatic  sea  clutter,  as  discussed  in  Section  2.2.5  ,  where  the 

“standard  atmosphere”  propagation  factor  is  to  be  removed  and  the  grazing  angles  are  to  be  calculated 
using  TEMPER,  would  not  be  straightforward.  Eor  these  reasons,  the  computational  steps  required  to 
handle  ducting  conditions  are  not  included  in  the  model. 

3  COMPARISON  WITH  MEASURED  DATA 

Comparisons  with  measured  data  are  essential  to  corroborate  the  theoretical  model.  Currently,  limited 
data  is  available  for  comparison  with  the  bistatic  model  especially  in  the  glistening  region.  The  most  well- 
known  are  the  experiments  carried  out  by  Beard  [15]  in  1961  between  two  antennas  of  equal  height  and 
about  3.5  km  apart.  Recently  Haspert  [16]  has  also  made  some  low  grazing  angle  measurements  and 
compared  his  measured  results  to  theory  discussed  by  Barton  [17].  Barton’s  theory  is  also  based  on 
Beckmann  and  Spizzichino.  The  biggest  issue  in  making  multipath  measurements  from  the  glistening 
zone  is  the  small  relative  delay  between  the  direct  and  multipath  returns.  These  short  delays  demand  use 
of  narrow  pulses  with  low  sidelobes  to  enable  resolution  between  direct  and  indirect  returns.  These 
narrow  pulses  need  large  bandwidths.  The  small  differences  in  velocities  between  direct  and  indirect 
returns  makes  using  Doppler  difficult. 

However,  there  is  some  measured  data  available  in  the  wide-angle  or  backscatter  region.  The  next 
section  discusses  comparisons  with  some  examples  of  this  data.  There  are  no  discussions  of  specular 
measurement  comparisons  as  the  developed  model  is  primarily  concerned  with  bistatic  reflectivity  results. 


3.1  Pidgeon  Data  Set 


In  a  1965  paper  by  V.W.  Pidgeon  [12],  an  experimental  program  was  carried  out,  where  the  bistatic 
radar  cross  section  was  measured  as  a  function  of  transmitter  and  receiver  depression  angles.  In  Eigure  14, 
two  different  lines  are  drawn.  The  “blue”  curve  (parametrically  calculated  by  the  authors)  is  the  drawn 
over  one  of  Pidgeon’s  data  and  its  equation  is: 

0-0=  -47  +  101ogio^5 


(23) 
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where  is  the  transmitter  depression  angle.  This  parametric  equation  is  valid  only  for  sea  state  3  and 
vertical  polarization.  The  NRL  histatic  model  is  the  basis  for  the  “red”  curve  for  increasing  transmitter 
depression  angle  at  sea  state  3,  vertical  polarization  and  at  C-Band. 

The  two  curves  show  good  agreement  between  the  measured  and  modeled  data.  For  example,  at 
equal  to  0.6°,  the  experimental  data  predicts  the  bistatic  RCS  to  be  about  -50dB  while  the  NRL  model 
predicts  this  value  to  be  about  -52dB. 
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Figure  14:  Summary  of  Bistatic  Cross-section  for  Various  Sea  and  Wind  conditions;  Transmitter 

Polarization:  vertical;  Frequency:  C-band 


3.2  Kochanski  Data  Set 

A  more  recent  paper  (1992)  by  T.  P.  Kochanski  [13],  at  M.I.T  Lincoln  Laboratory,  summarizes 
measurements  of  bistatic  sea  backscatter  at  X-band  (lOGHz)  using  the  Airborne  Seeker  Test  bed  and  a 
CW  (continuous  wave)  surface  based  illuminator.  These  measurements  were  made  at  sea  state  one  and  at 
vertical  polarization.  The  transmitting  antenna  had  a  fixed  depression  angle  while  the  depression  angle  of 
the  receiving  antenna  located  on  an  aircraft  varied.  The  transmitter  depression  angle  remained  fixed  at 
0.3°.  The  receiver  depression  varied  from  5°  to  40°.  The  experimental  data  showed  that  the  bistatic  RCS 
was  essentially  independent  of  receiver  depression  angle. 

Figure  15  shows  the  Kochanski  experimental  data  overlaid  with  the  results  obtained  using  the  NRL 
bistatic  model.  We  observe  the  same  independence  from  the  transmitter  depression  angle  in  the  NRL 
bistatic  model.  From  the  results,  the  average  experimental  bistatic  RCS  is  about  -50dB  while,  the  NRL 
bistatic  model  puts  this  value  at  about  -58dB.  Once  again,  discounting  for  experimental  errors,  there  is  a 
reasonable  match  between  experimental  and  modeled  data. 
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Depression  Angle  (deg) 

Figure  15:  Comparing  Kochanski’s  Experimental  Results  at  X-band  with  Results  Obtained  using  NRL 

Bistatic  Model 

3.3  Ewell  Data  Set 

This  1982  publication  by  G.  W.  Ewell  [14]  describes  experiments  carried  out  at  X-band  for  both 
horizontal  and  vertical  polarization  for  a  range  of  transmitter  and  receiver  depression  angles.  The  results 
plots  the  bistatic  and  monostatic  RCS  vs.  the  bistatic  angle,  defined  as  the  angle  between  the  receiver,  the 
transmitter  and  the  “aim  point”  between  the  two  antennas  as  shown  in  figure  16.  The  data  presents  the 
ratio  of  the  bistatic  a°  values  and  the  monostatic  a°  values.  Absolute  values  were  not  given,  maybe 
because  of  an  unreliable  radar  system  calibration. 

To  get  the  data  shown  in  figure  17,  the  transmitting  antenna  has  a  look  direction  of  90°  in  sea  state  4. 
The  separation  distance  between  the  receiving  and  transmitting  antenna  is  at  3.5km.  The  height  of  the 
transmitting  antenna  is  22.9m  while  the  height  of  the  receiving  antenna  is  at  8.2  m.  Horizontal  and 
vertical  data  was  then  taken  for  bistatic  angles  ranging  from  about  20°  to  60°.  Figure  17  shows  the  data 
calculated  using  the  NRL  bistatic  model  overlaid  on  the  Ewell  results  for  both  the  horizontal  and  vertical 
polarizations.  For  the  three  cases  shown,  the  bistatic  RCS  is  smaller  than  the  monostatic  RCS  by  about  10 
dB.  This  is  as  expected  since  the  monostatic  RCS  is  based  on  the  transmitter  depression  angle  while  the 
bistatic  RCS  is  based  on  the  smaller  of  the  receiver  and  transmitter  depression  angle. 

However,  the  results  show  that  for  most  bistatic  angles,  the  ratio  of  the  bistatic  RCS  to  the  monostatic 
RCS  from  the  NRL  model  never  falls  below  -20dB,  which  is  greater  than  the  experimental  data  especially 
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at  the  larger  bistatic  angles.  In  the  model,  for  the  geometry  shown  in  Fig.  6,  Eq.  (21)  can  be  used  to 
calculate  both  the  monostatic  and  bistatic  RCS  with  different  a  values.  For  the  monostatic  case,  a  is  equal 
to  the  transmitter-grazing  angle  while  for  the  bistatic  case,  a  is  equal  to  the  smaller  of  the  transmitter  or 
the  receiver-grazing  angle.  For  most  scenarios,  the  receiver  is  the  smaller  grazing  angle.  However,  for  no 
geometry  is  the  receiver  grazing  angle  so  much  smaller  than  the  transmitter-grazing  angle  such  the  ratio 
of  the  bistatic  to  the  monostatic  RCS  falls  much  less  than  10  dB  for  fixed  frequency,  polarization  and  sea 
state.  Hence,  further  analysis  may  be  necessary  to  understand  the  discrepancy  between  the  measured  and 
modeled  data  for  this  case. 


X 


Figure  16:  Geometry  and  Nomenclature  for  the  Ewell  Experiment  [14] 
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Figure  17:  Comparison  of  Ewell  Experimental  Data  with  the  NRL  Bistatic  Model  Data 
4  CONCLUSIONS 


Modeling  of  bistatic  sea  surface  reflectivity  is  a  complex  problem  due  to  the  many  parameters 
defining  the  geometry  of  interest.  The  current  state  of  the  art  invokes  three  different  components:  coherent 
scatter  (reflection)  from  the  specular  region,  strong  directive  diffuse  scatter  from  the  glistening  region, 
and  wide  angle  diffuse  scatter  in  directions  outside  the  glistening  region.  The  report  describes  in  detail  the 
steps  used  to  model  the  diffuse  scatter  assuming  a  Gaussian  sea  surface  model.  This  basis  of  this 
geometrical  model  is  primarily  the  work  done  by  Beckmann  and  Spizzichino  in  the  1960s.  Currently, 
there  is  a  serious  lack  of  a  mature  bistatic  sea  reflectivity  model  because  of  limited  measured  bistatic  sea 
clutter  data  to  seriously  anchor  any  modeling  attempts.  There  are  still  many  unknown  issues  related  to 
monostatic  sea  clutter,  for  example  sea  spikes,  which  will  certainly  impact  bistatic  clutter  modeling,  and 
which  are  still  very  much  under  active  research.  The  model  presented  here  represents  the  best  available 
theory  but  uncertainties  must  nevertheless  be  associated  with  any  quantitative  results  obtained. 
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Appendix  A  -  MATLAB  IMPLEMENTATION  OE  BISTATIC  MODEL 

A  MATLAB  function  was  written  to  compute  the  bistatic  sea  reflectivity  for  a  given  set  of 

parameters.  This  function  will  calculate  the  average  horizontal  and  vertical  bistatic  reflectivity  for  a 
specified  surface  patch.  With  reference  to  the  program  listing  shown  below  as  Figure  1  the  inputs  are  as 
follows: 

1.  Height  of  the  receiving  antenna  (/i^)  in  meters 

2.  Angle  of  elevation  (0£j)  to  the  transmitting  antenna 

3.  Sea  State  (SS) 

4.  Distance  between  the  receiving  and  transmitting  antennas  (D)  in  meters 

5.  Frequency  of  operation  (FGHz)  in  GHz 

6.  The  x-coordinate  of  the  surface  patch  in  meters 

7.  The  y-coordinate  of  the  surface  patch  in  meters 

8.  Indicate  whether  shadowing  is  to  be  applied  or  not 

9.  Indicate  the  polarization  of  the  incident  ray  (choose  between  Horizontal  or  Vertical) 

The  outputs  of  this  code  are  the  following: 

1.  The  average  co-pol  and  cross-pol  surface  reflectivity  for  the  given  x-  and  y-coordinates 
in  dB 

2.  The  angle  from  the  receiver  horizontal  to  the  surface  patch 

3.  The  angle  from  the  transmitter  horizontal  to  the  surface  path  ((prx) 

4.  The  off-axis  angle  from  the  receiver  to  the  surface  patch  (0i) 

5.  The  off-axis  angle  from  the  transmitter  to  the  surface  patch  (82) 

6.  Grazing  angle  from  receiver  (y^) 

7.  Grazing  angle  from  transmitter  (yj-) 

8.  Length  of  path  from  receiver  to  surface  patch  (R^)  in  meters 

9.  Length  of  path  from  transmitter  to  surface  patch  (R2)  in  meters 

10.  Direct  path  length  between  receiver  and  transmitter  (R^)  in  meters 

11.  Height  of  transmitting  antenna  (hj-)  in  meters 

The  main  report  defines  the  equations  used  in  this  MATLAB®  listing. 
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function[sigmaCoPol_dB,  slgmaXPol_dB, . . . 

grazRx,  grazTx,  phiRx,  phiTx,  thetaRx,  thetaTx,  Rl,  R2,  Rd,  hT]  =  ... 

ReflectivityCoeff_Calculation(hR,  thetad,  SeaState,  D,  FGHz,  xPatch,  yPatch,  Shadowing,  TxPol,  Type,hT) 
cic; 

%%  CONSTANTS 
Deg2Rad  =  pi/180; 

SpeedofLight  =  3e8; 

Re  =  8500e3; 

if  nargin  ==  0 
hR  =  20; 

D  =  10e3; 
thetad  =  3.0; 

SeaState  =  3; 

xPatch  =  50; 
y  Patch  =  0; 

FGHz  =  3; 

Shadowing  = 

TxPoi  =  V'; 

Type  =  1; 

end; 

%Define  siopefrom  Sea  State 

if  SeaState  ==  0 
tanbetaO  =  0.05; 
eiseif  SeaState  ==  1 
tanbetaO  =  0. 12; 
eiseif  SeaState  ==  2 
tanbetaO  =  0. 14; 
eiseif  SeaState  ==  3 
tanbetaO  =  0. 15; 
eiseif  SeaState  ==  4 
tanbetaO  =  0.16; 
eiseif  SeaState  ==  5 
tanbetaO  =  0. 13; 
eiseif  SeaState  ==  5 
tanbetaO  =  0.22; 
eiseif  SeaState  ==  7 
tanbetaO  =  0.25; 
end; 

%Eiectricai  properties  of  sea  water 
iambda  =  SpeedofLight/(FGHz*le9); 
epsrc  =  80-  60  *  sgrt(-l)  *  iambda  *  4; 
mure  =  1; 
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y Earth  =  sqrt(epsrc/murc); 


%%  START  CALCULATION  OF  SURFACE  REFLECTIVITY 
alpha  =  D/Re; 

TD  =  thetad  *  Deg2Rad; 

%Calculate  the  height  of  Transmitter 
if  Type  ==  1 

beta  =  pi  -  (TD  +  pi/2)  -  alpha; 

hT  =  ((Re  +  hR)  *  sin(TD  +  pi/2)/ (sin(beta)))  -  Re; 

end; 

Rd  =  sqrt((hT  -  hR)''2  +  4*(Re+hR)*(Re+hT)*sin(alpha/2)''2); 

xl  =  sqrt(xPatch''2  +  yPatch''2); 
alphal  =  xl/Re; 

ifxPatch  ==  0 
theta4  =  0; 
elseifxPatch  <  0 
theta4  =  pi/2  +  acos(xPatch/xl); 
else 

theta4  =  acos(xPatch/xl); 
end; 

thetaRx  =  theta4/Deg2Rad; 

x2  =  sqrt(xl  ''2  +  D''2  -  2  *  xl  *  D  *  cos(theta4)); 
thetas  =  acos((x2''2  +  D''2-xl ''2)/(2  *x2*  D)); 

thetaTx=  theta3/Deg2Rad; 


alpha2  =  x2/Re; 

R1  =  sqrt(hR''2  +  4  *  Re  *  (Re  +  hR)  *  sin(alphal/2)''2); 

R2  =  sqrt(hT''2  +  4  *  Re  *  (Re  +  hT)  *  sin(alpha2/2)''2); 

%Calculate  grazing  and  incidence  angles  from  receiver  hR 
verthl  =  (hR  +  Re)  -  (Re/cos(alphal)); 
horzhl  =  Re  *  tan(alphal); 

grazl  =  acos((Rl''2  +  horzhl''2  -  verthl''2) / (2  *  R1  *  horzhl)); 


if  isnan(grazl) 
grazl  =  pi/2; 
end; 

%Grazing  angle  from  reflecting  surface  to  Receiver 
grazRx  =  grazl/Deg2Rad; 

thetal  =  pi/2  -  grazl; 
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%Angle  from  Rx  horizontal  to  Reflecting  surface 

phiRx  =  90  -  acos(((Re+hR)''2  +  Rl'^l  -  Re''2)/(2  *  R1  *  (Re+hR)))/Deg2Rad; 

verth2  =  (hT  +  Re)  -  Re/(cos(alpha2)); 
horzh2  =  Re  *  tan(alpha2); 

graz2  =  acos((R2''2  +  horzh2''2  -  verth2''2)/(2  *  R2  *  horzh2)j; 

if  isnan(graz2) 
graz2  =  pi/2; 
end; 

%Grazing  angle  from  reflecting  surface  to  Transmitter 
grazTx  =  graz2/Deg2Rad; 

theta2  =  pi/2  -  graz2; 

%Angle  from  Tx  horizontal  to  Reflecting  surface 

phiTx  =  90  -  acos(((Re  +  hTj''2  +  R2''2  -  Re''2)/(2  *  R2  *  (Re  +  hT)))/Deg2Rad; 


if  y Patch  <0 

totAngle  =  abs(theta4)  +  theta3; 
else 

totAngle  =  2  *  pi  -  (abs(theta4)  +  thetaS); 
end; 

tanbeta  =  sgrt(sin(thetal)''2  -  2  *  sin(thetal)  *  sin(theta2)  *  cos(totAngle)  +  sin(theta2)''2) / ... 
(cos(thetal)  +  cos(theta2)); 

if  strcmpi(Shadowing,'Y')  ==1 
vl  =  abs(cot(thetal))/(sgrt(2)*tanbeta0); 
m2  =  abs(cot(tbeta2))/(sgrt(2)*tanbeta0); 

B1  =  (exp(-vl  *vl)./vl-  sgrt(pi)  *  erfc(vl)) ./  (2  *  sgrt(pi)); 

B2  =  (exp(-M2  *  m2)  ./  m2  -  sgrt(pi)  *  erfc(M2)) ./  (2  *  sgrt(pi)); 

Shadow_Factor  =  (1  +  erf(Ml))  *(1+  erf(M2))  ./(4.*(1+B1+  B2)); 
else 

Shadow_Factor  =  1; 
end; 

if  Mertbl  >  0  &&  Mertb2  >  0 

tempi  =  (l+tanbeta''2)''2/(tanbeta0''2)  *  exp(-(tanbeta/tanbeta0)''2); 
tempi  =  tempi  *  Shadow_Factor; 

%Calculate  alpha  -  to  include  Polarization 

alpha  =  acos((sgrt(l  -  cos(grazl)*cos(graz2)*cos(totAngle)  +  sin(grazl)*sin(graz2)))/sgrt(2)); 
Alpha  =  pi/2  -  alpha; 

ifstrcmpi(TxPol,  'H')  ==1 

RHoriz  =  (sin(Alpha)  -  sgrt(YEarth^2  -  cos(Alpha)''2)) ... 
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/(sin(Alpha)  +  sqrt(YEarth''2  -  cos(Mpha}''2}}; 
elseif  strcmpi(TxPol,  V'j  ==1 

RVert  =  (YEarth''2  *  sin(Alpha)  -  sqrt(YEarth''2  -  cos(Alpha)''2)) ... 

/(YEarth''2  *  sin(Alpha)  +  sqrt(YEarth''2  -  cos(Alpha)''2jj; 

end; 

cosX  =  cos(thetal)  *  cos(theta2)  -  sin(thetal)  *  sin(theta2)  *  cos(totAngle); 
sinX  =  sqrt(l-cosX''2); 

sinbetal  =  sin(thetal)  *  sin(totAngle)/sinX; 
sinbeta2  =  sin(theta2)  *  sin(totAngle)  / sinX; 

cosbetal  =  (sin (theta 2}  *  cos(thetal)  +  cos(theta2)  *  sin(thetal)  *  cos(totAngle))/sinX; 
cosbeta2  =  (sin(thetal)  *  cos(theta2)  +  cos(thetal)  *  sin(theta2)  *  cos(totAngle))/sinX; 

ifstrcmpi(TxPol,  'H'j  ==  1 

sigmaCoPol  =templ  *  real((abs(RHoriz)  *  cosbetal  *  cosbeta2)''2); 
sigmaXPol  =  tempi  *  real((abs(RHoriz)  *  cosbetal  *  sinbeta2)''2); 

elseif  strcmpi(TxPol,  V'j  ==  1 

sigmaCoPoi  =  tempi  *  real((abs(R\/ert)  *  cosbetal  *  cosbeta2)''2); 
sigmaXPol  =  tempi  *  real((abs(R\/ert)  *  cosbetal  *  sinbeta2)''2); 
end; 


%Have  to  include  Monostatic  for  wide  angle  scattering 
%  ifxPatch  <  0 
%  graz  =  grazRx; 

%  else 

graz  =  min(grazRx,grazTx);  %Choose  the  smaller  of  the  two  grazing  angles 
%  end; 

ifstrcmpi(TxPol,  'H')  ~1 

CCl  =  -73.0;  CC2  =  20.781;  CC3=  7.351;  CC4=25.65;  CC5  =  0.0054; 

elseif  strcmpi(TxPol,  V)  ~1 

CCl  =  -50. 796;  CC2  =  25.93;  CC3  =  0.7093;  CC4  =  21.588;  CC5  =  0.00211; 

end; 

temp2  =  CCl  +  CC2  *  loglO(sin(graz*Deg2Rad))  +  ... 

(27.5  +  CC3  *  graz)  *  logl0(FGHz)/(1.0  +0.95  *  graz)  +  ... 

CC4  *  (SeaState  +  1) ''  (l/(2.0  +  0.085  *  graz  +  0.033  *  SeaState))  +  CC5*graz.''2; 

temp2  =  10''(temp2/10); 

A  =  [sigmaCoPol,  temp2]; 

sigmaCoPol  =  max(A); 

A  =  [sigmaXPol,  0]; 

sigmaXPol  =  max(A); 
else 

SigmaCoPol  =  0; 

SigmaXPol  =  0; 
end; 
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%Co-  and  X-Pol  Bistatic  RCS  in  dB 

sigmaCoPol_dB  =  10.  *loglO(sigmaCoPol); 
sigmaXPol_dB  =  10.  *loglO(sigmaXPol); 

ifsigmaCoPol_dB  <=  -80 
sigmaCoPol_dB  =  -80; 
end; 

ifsigmaXPol_dB  <=  -80 
sigmaXPol_dB  =  -80; 
end; 


Figure  A.l:  MATLAB  Listing  for  Calculation  of  Average  Sea  Reflectivity  using  Geometric  Bistatic  Model 
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Appendix  B  -  MATLAB  GRAPHICAL  USER  INTERFACE  (GUI)  FOR  BISTATIC 
CALCULATIONS 

The  MATLAB  function  described  in  Appendix  A  was  used  to  construct  a  Graphical  User  Interface 
(GUI)  to  provide  a  convenient  way  to  use  the  code  for  a  given  surface  geometry.  This  MATLAB  program 
determines  the  basic  reflectivity,  the  bistatic  RCS,  as  well  as  received  power  and  time  delay  distribution. 

The  inputs  to  the  GUI  include: 

•  Frequency  of  operation  in  GHz 

•  Transmit  and  receive  polarization 

•  Range  in  km  between  receiver  and  transmitter 

•  Height  of  receiver  in  meters 

•  Elevation  angle  to  transmitter 

•  Gain  and  sidelobe  level  of  receiving  antenna,  which  has  a  Gaussian  shape  and  a  constant  sidelobe 
level.  The  transmit  antenna  is  isotropic 

•  Sea-state  can  be  chosen  as  shown  in  Table  1 

•  The  calculations  are  based  on  the  user  input  of  flat  or  spherical  earth 

•  The  option  of  including  shadowing  and/or  wide-angle  scatter  is  available 

•  The  resolution  of  the  plot  as  well  the  size  of  the  surface  can  also  be  specified 

The  code  outputs  the  maximum  scattered  power  and  determines  the  diffused  scattering  coefficient.  It 
also  plots  as  output 

•  The  glistening  surface 

•  The  average  bistatic  RCS  depending  on  unit  patch  size 

•  The  apparent  bistatic  RCS  which  includes  the  effect  of  gain  of  the  main  beam  of  the  receiving 
antenna, 

•  The  average  diffused  power  from  the  surface  depending  on  individual  patch  sized  defined  by  user 

•  Time  vs.  power  plots  which  indicates  the  delay  of  the  diffused  power  from  the  surface  as  it 
reaches  the  receiving  antenna.  This  is  very  important  when  determining  the  degradation  in  the 
sidelobe  canceller  because  scattered  power  with  high  delays  cannot  be  easily  cancelled  by  an 
adaptive  canceller 

•  Range  vs.  power  plots,  which  allow  the  user  to  determine  at  which  range  most  diffused  power 
comes  from. 
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Figure  B.l  shows  the  GUI  for  the  MATLAB®  code  and  Figure  B.2  shows  examples  of  results 
produced  hy  the  code.  All  code  is  available  from  NRL  upon  request. 


Figure  B.l:  MATLAB  GUI  (Graphical  User  Interface)  to  determine  diffuse  power  from  a  planar  surface  for  a 

given  geometry 
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X-direction  (km)  X-dIrection  (km) 


•  This  plot  shows  the  histatic  RCS  over 
the  entire  surface  of  interest 

•  The  equation  used  for  this  plot  is 
Equation  (10)  of  the  main  report 

•  The  plot  is  limited  due  to  the  horizon 


of  Surface 

In  this  plot,  the  gain  of  the  transmit 
and  receive  antennas  are  also  included 

Most  of  the  diffused  power  comes 
from  close  to  the  receive  antenna 
where  the  delays  can  he  quite  long 


Time  (ns) 


(b)  Diffuse  Power  over  Given  Surface 

•  This  is  cumulative  distribution  plot  of 
the  incoming  power  vs.  time  (ns) 

•  The  power  numbers  have  been 
normalized  to  the  maximum  incoming 
power 

•  In  this  case  more  than  50%  of  the 
power  has  delays  greater  than  3ns 
(since  no  wide-angle  scatter)  has  been 

(c)  Plot  of fractional  power  vs.  time  (ns) 
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9^  =  1.2  and  Jammer  Azimuth  Angles  o 

•  This  is  cumulative  distribution  plot 
of  the  incoming  power  vs.  range 
(km) 

•  The  power  numbers  have  been 
normalized  to  the  maximum 
incoming  power 

•  This  plot  allows  the  user  to 
understand  the  range  from  which 
most  of  the  power  is  coming  from 

Range  (km)  ^  ® 

(d)  Plot  of  fractional  power  vs.  range  (km) 

The  tab  gives  a  summary  of  the  results 
calculated  by  the  program 

Specular  coefficient 
Diffuse  scattering  coefficient 

Maximum  multipath  power 
(e)  Summary  of  results 

Figure  B.2:  Example  of  results  available  from  “Bistatic  Sea  Surface  Reflectivity”  Code 


SUMMARY  RESULTS. 


Specular  ScaRenng 

Diffuse  Scattenng 

CoeJ  0.284 

Coef 

0.00142 

Multipath  Power  in  Mam  Beam 

(06)  "M 

(Relsove  to  Direct  Petti  end  IsotrofiK  Receive  Antenne) 
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